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O ■ 

o : 

^ ' We investigate the effect that the structure of GRB jets has on the afterglow 

Q ■ hght curves for observers located at different viewing angles, 6'obs, from the jet 

Q ■ symmetry axis. The largest uncertainty in the jet dynamics is the degree of 

l/^ ■ lateral energy transfer. Thus, we use two simple models, that make opposite and 

^ ■ extreme assumptions for this point, and calculate the light curves for an external 

^ ■ density that is either homogeneous, or decreases as the square of the distance from 

Q ' the source. The Lorentz factor and kinetic energy per unit solid angle are initially 

^ ■ taken to be power laws of the angle from the jet axis. We perform a qualitative 

^ ■ comparison between the resulting light curves and afterglow observations. This 

^ ■ constrains the jet structure, and poses problems for a 'universal' jet model, where 

O ■ all GRB jets are assumed to be intrinsically identical, and differ only by our 

^ ■ viewing angle, 6'obs. 

2 ' Subject headings: gamma rays: bursts — gamma-rays: theory — shock waves — 

^ ■ relativity — radiation mechanisms: nonthermal 

> ' 

' 1. Introduction 

There are several different lines of evidence in favor of collimated outflows, also referred 
to as jets, in gamma-ray bursts (GRBs). Perhaps, the best evidence so far is the achromatic 
break seen in the light curves of many (though not all) GRB afterglows. A collimated outflow 
also helps to reduce the estimate for the total energy output in gamma-rays that is inferred 
from the fiuence in GRBs with known redshifts, which in some cases approaches and in 
one case (GRB 990123) even exceeds the rest energy of a solar mass star, for a spherically 
symmetric emission. Such an energy output in gamma-rays is hard to produce in any model 
that involves a stellar mass object. Furthermore, a non-spherical flow can also manifest itself 
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through hnear polarization (Sari 1999; GhiseUini & Lazzati 1999) as was indeed observed for 
a few afterglows (Covino et al. 1999; Wijers et al. 1999; Rol et al. 2000). 

Most GRB jet models consider an outflow that is uniform within some finite, well 
defined, opening angle around its symmetry axis, and where the Lorentz factor, energy 
density etc. drop sharply beyond this opening angle (Rhoads 1997, 1999; Panaitescu & 

Meszaros 1999; Sari, Piran & Halpcrn 1999; Kumar & Panaitescu 2000; Moderski, Sikora & 
Bulik 2000; Granot et al. 2001, 2002). Such a uniform jet with a sharp, well defined, edge 
shall be referred to as a 'top hat' jet. The possibility that GRB jets can display an angular 
structure, where the kinetic energy per unit solid angle, e, and the Lorentz factor, F, of the 
GRB outflow vary smoothly as power laws in the angle, 9, from the jet axis, was proposed 
by Meszaros, Rees & Wijers (1998). We shall refer to such smoothly varying relativistic 
outflows as 'structured' jets (as opposed to a 'top hat' jet, that has no inner structure). 

Recently, several different groups have analyzed afterglow observations within the frame 
work of the 'top hat' jet model, and have inferred a relatively narrow distribution both for 
the total energy output in gamma-rays (Frail et al. 2001) and in the initial kinetic energy 
of the relativistic outflow (Panaitescu & Kumar 2001; Piran et al. 2001). These results may 
alternatively be interpreted as GRB jets having a universal structure, which is intrinsically 
the same for all GRBs, and the observed differences between different GRBs are a result of 
different viewing angle, ^obs, w.r.t. the jet symmetry axis (Lipunov, Postnov & Prokhorov 
2001; Rossi, Lazzati & Rccs 2002; Zhang & Meszaros 2002). Whereas in the 'top hat' jet 
interpretation, the jet break time, tj, depends mainly on the initial opening angle of the jet, 
9j (and also has a smaller dependence on its energy per unit solid angle and on the external 
density, e.g. Sari, Piran & Halpern 1999), in the universal 'structured' jet interpretation, tj 
depends mainly on the viewing angle, ^obs; and the light curve is similar to that for a 'top 
hat' jet with an opening angle 9j — ^obs a-nd the same value of the energy per unit solid angle 
of the 'structured' jet at ^ = ^obs- 

While the evolution of 'top hat' jets and their hght curves have been widely investigated 
(Rhoads 1999; Panaitescu & Meszaros 1999; Kumar & Panaitescu 2000; Moderski, Sikora 
& Bulik 2000; Granot et al. 2002), including numerical simulations of the jet dynamics and 
calculation of the resulting light curves (Granot et al. 2001), much less work has been done 
on 'structured' jets. In an accompanying paper (Kumar & Granot 2002) we calculate the 
dynamics of 'structured' relativistic jets by solving relativistic fluid dynamics equations, and 
demonstrate that a simple analytic or semi-analytic model can qualitatively reproduce the 
results for the jet dynamics and for the light curves as long as the Lorentz factor along the 
jet axis is of order 4 or larger. In this paper we use two simple models for the jet dynamics, 
which are designed to bracket the true hydrodynamic evolution of the jet, to calculate the 
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afterglow light curves for a wide range of parameters for the structured jet and the external 
density. A qualitative comparison of the light curves with afterglow observations provides 
constraints on the jet structure and density profile in its immediate vicinity. 

In §2 we present the physical model. We begin with the initial conditions (§2.1), then 
describe our two simple dynamical models (§2.2) and finally outline the procedure for calcu- 
lating the light curves (§2.3). Our results are presented in §3, and in §4 we discuss our main 
conclusions. 



2. The Physical model 

The jet is assumed to posses axial symmetry, so that all the hydrodynamic quantities 
at a given lab frame time, t, depend only on the angle, 9, from the jet symmetry axis. The 
radial profile of the outflow is ignored in this simple treatment, and the shocked material is 
approximated by a thin shell that is located at the same radius, R, of the shock. 

The observed afterglow fight curves depend on the angular structure of the GRB jet, 
its hydrodynamic evolution, and the viewing angle, ^obs, w.r.t. the jet symmetry axis. One 
may therefore use the shape of the afterglow light curves to constrain the angular structure 
of the jet, as well as to infer the viewing angle, ^^obs- For a 'top hat' jet of opening angle 
^obs < ^0 is required in order to observe the prompt gamma-ray emission. In this case the 
differences in the light curves between different 6'obs in the range < 6'obs < are rather 
small (Granot et al. 2001, 2002). However, for structured jets we expect large differences in 
lightcurves for different Ogbs- 

2.1. The Initial Conditions 

The initial conditions are chosen at a lab frame time, to, for which on the one hand, the 
internal shocks have ended, and on the other hand, no significant deceleration due to the 
sweeping up of the external medium has yet occurred (i.e. Ris{0) < R{9,to) < -Rdec(^) for 
all relevant 9). For simplicity, the original ejecta is assumed to remain cold, even near the 
deceleration radius, where it is expected to be heated by the passage of the reverse shock. 
This approximation might introduce inaccuracies of order unity near the deceleration epoch, 
corresponding to an observed time < a few minutes for a typical GRB, but it should have 
no effect on lightcurves at later times. Energy conservation implies 

e{9, to) = i^o{0, to) [T{9, to) - 1] + i^s{9, to) [T\9, to) - l] ^ fMo{0, to)T{9, to) , (1) 
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where hq and fig are the rest mass per unit solid angle of the original ejecta and of the swept 
up external medium, respectively, e is the energy (not including the energy associated with 
the rest mass) per unit solid angle in the outflow, F is the bulk Lorentz factor of the shocked 
material, and to is chosen so that iJ,s{d,to)T{9,to) <^ /io(^, ^o) for all 9. 

The initial kinetic energy per unit solid angle, and initial Lorentz factor (minus 1) are 
assumed to be power laws in 9, outside of a core of opening angle 9^- 

e(^,to) = eoe-« , r(^,to) = i + (ro-i)G-\ (2) 

where cq and Fq are the initial kinetic energy per unit solid angle and Lorentz factor at the 
jet axis, and 



The velocity is assumed to be in the radial direction, so that the initial radius is given by 



R{9, to) = ioVl-r-2(^,io) , (4) 

and the lateral transfer of matter can be neglected^. The external mass density profile is 
assumed to be a power law in the distance r from the source, Pextij) = Ar^^ . This implies 



f^^^'^^ , AR{9,tf-^ 
f,o{9, t) = f,o{9, to) = " - [ri9, to) + 1] i^siO, to) . (6) 

U - -Lj 



2.2. The Jet Dynamics 

The main uncertainty in the jet dynamics is the degree of lateral transport of energy, 
from small to large angles, 9, w.r.t. the jet axis. We therefore make two alternative and 
extreme assumptions regarding this transport: in model 1 we assume that the energy per 
unit solid angle, e, does not change with time, e{9, t) = t{9, to), while in model 2 we assume 
the maximal averaging of e over the angle 9 that is consistent with causality. This is done 
by averaging over the initial distribution of e. 



more rigorous treatment of the jet hydrodynamics (Kumar & Granot 2002) shows that these assump- 
tions are reasonable. 
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where 9- (^+) is the angle below (above) 6 (i.e. 9- < 9 < 9+), to which a hypothetical 
sound wave that originated at 9 at the initial time to would have propagated to. Initially, 
^-(^o) — ^+(^o) — ^, and these angles subsequently evolve according to 



d9i 
'dt 



= ±- 



Cs{9±,t) 



7±, 



t)R{9±,t) 



/ 7(7-i)[r-i] 
i + 7[r-i] 



(8) 



where is the local sound speed, and 7 — {AT + l)/3r is the adiabatic index. This simple 
local scheme does not, in general, conserve the total energy in the outflow. The global energy 
conservation is imposed by renormalizing e. 



£^^d9'sm9'e{9',t) 



(9) 



For both models 1 and 2, the radius of the thin shell of matter changes as^ 

dR{9,t) 



dt 



(10) 



while lis and /iq are given by equation (5) and (6), respectively, and the Lorentz factor of 
the shocked material, F, is obtained by solving the equation for energy conservation. 



e = (F - l)/io + (r' - l)fis 



JM)_ 



^ ^ Ansje + Ho + Us) _ ^ 



(11) 



2.3. The Light Curves 

The local emissivity is calculated using the conventional assumptions of synchrotron 
emission from relativistic electrons that are accelerated behind the shock into a power law 
distribution of energies, A^(7) oc 7"^ for 7 > where the electrons and the magnetic field 
hold fractions eg and e^, respectively, of the internal energy. The shape of the local spectral 
emissivity is approximated as a broken power law with breaks at the typical synchrotron 
frequency Um and at the cooling frequency Uc (Sari, Piran & Narayan 1998). As our main 
focus in this work is the effect of the jet dynamics on the afterglow emission, we keep the 
expression for the local emissivity fairly simple, and leave the inclusion of additional features 
such as the self absorption frequency, inverse Compton emission and its effect on the electron 



^There is a small difference between the velocity of the shock front and that of the fluid just behind the 
shock (Blandford & McKee 1976), however we neglect this distinction in the present treatment, for simplicity. 
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cooling, etc., for later applications, since such complications might make it hard to pinpoint 
the effects of the jet dynamics on the light curves. 

The hght curves for observers located at different angles, ^obsi w.r.t. the jet axis, are 
calculated by applying the appropriate relativistic transformation of the radiation field from 
the local rest frame of the emitting fluid to the observer frames, and integrating over equal 
photon arrival time surfaces (Granot, Piran & Sari 1999; Kumar & Panaitescu 2000; Granot 
et al. 2001). The radiation calculation is in this sense rigorous and takes into account all 
relevant effects, so that the resulting light curves accurately reflect what is expected for a 
given jet structure and dynamics, where the latter are much less certain. 



3. Results 

In Figures 1-3 we show the evolution of the hydrodynamic quantities according to the 
two models described in the previous section, for a constant density external medium (A; = 0). 
In Figure 1 we show the evolution of the Lorentz factor for (a, b) = (0, 2). For a = 0, models 1 
and 2 become identical, so that this flgure applies to both models. Since e{6, t) = = const, 
the evolution after the deceleration time is the same for all 6, and different 6 differ just by 
the deceleration time, tdec, and the corresponding deceleration radius, Rdec, that are given 
by Us ~ /^o/r ~ e/r^, or 



-^dec 



(3 - k)eo 

Ac^rl 



(3-fe) 2b- g i?.„„ 2!,(4-fc)-a 

, tdec = t . « e (3-*) . 12 

' 2cr2(io) ^ ^ 



As can be seen from equation (12), for a — 0, the deceleration occurs flrst at small 9, and 
gradually proceeds to larger angles. This can be nicely seen in Figure 1. 

This situation is similar for (a, 6) = (2,2), but is reversed for (a, 6) = (2,0), as can be 
seen in Figures 2 and 3, where in the latter case the deceleration occurs flrst at large angles 
9 and then proceeds to smaller angles. This is in agreement with equation (12). The upper 
panel of Figures 2 and 3 shows the energy per unit solid angle, e, which is initially the same 
for the two models, and remains unchanged for model 1, while for model 2 e decreases at 
small angles 9, and increases at large angles. This indicates a lateral transfer of energy from 
small to large angles. At late times, as the flow becomes sub-relativistic, e becomes almost 
independent of angle 9, and the flow approaches spherical symmetry, as is indeed expected 
to occur physically, since the spherically symmetric Sedov- Taylor self similar solution should 
be asymptotically approached in the non-relativistic regime. The middle and bottom panels 
show the Lorentz factor minus one, F — 1, for models 1 and 2, respectively. For model 1, 
it can be seen that while the flow is still relativistic. The Lorentz factor soon settles into 
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a profile of F — 1 oc 9"^ a,t 9 > 9c, instead of the initial F — 1 oc 9~'^ or F = const. This 
occurs since the fraction of the energy in the original ejecta soon becomes negligible implying 
e ~ /^s(r^ — 1) so that as long as the flow is relativistic, /i^ oc R^~^ k, const, since Rk, ct for 
all 9, and F — 1 oc e^/^ oc 9^"^. When the flow becomes non-relativistic, V"^ — \ ^ 2(F — 1) ^ jS"^ 
so that F — 1 oc 9^'^^^^^^\ which is still quite close to for k between and 2. For model 
2, a similar effect is seen at early times, when e{9,t) ~ e{9,to), but it becomes much more 
homogeneous at later times, as e{9) becomes more uniform. 

The light curves of models 1 and 2 for (a, 6) = (0,0), (0,2), (2,0), (2,2) and a uniform 
external medium {k = 0) are shown in Figure 4, while the temporal decay slope, a (defined 
by -Fj, oc t°), for the same light curves, is shown in Figure 5. We have added the spherical 
case, (a, 6) = (0,0), for comparison. For (a, 6) = (0,2), the light curve initially rises, before 
the deceleration time, tdec, which occurs at latter times for larger viewing angles ^obs, in 
accord with equation (12), when keeping in mind that the emission is dominated by ^ ~ ^obs 
at early times due to the relativistic beaming of the radiation emitted from the jet. After a 
time t, deceleration has occurred at 6* < 9dec(t), where 9dec is given by idec[^dec(^)] = t, and 
the light curves for ^^obs < ^dec(^) approach the spherical case, and hence become very close 
to one other. Keeping a = 0, e and Fq fixed, and increasing b from zero, tdcc(^obs) begins 
to shift to larger times, and the ratio of the deceleration times for two given values of ^obs 
grows as b increases. But after the time idec(6'obs), the hght curves still approach the same 
spherical hght curve, for all values of b. There is no jet break in the hght curve for a — 0. 
Since jet breaks are observed in many afterglows, this type of jet structure can be ruled out 
as a universal model for GRB jets. 

For (a, b) = (2, 0) and (2, 2) we find a jet break in the light curve at roughly the same 
time as predicted in previous works (Rossi, Lazzati & Rees 2002; Zhang & Meszaros 2002), 
i.e. when F(^ = ^obs) ~ ^^s- However, we also find some new and interesting features. For 
(a, b) — (2, 2) with Fq < 10^, the initial Lorentz factor at large angles, 9, is quite modest, 
resulting in a large deceleration time, tdec which can be as large as tdec ^ 1 day for 6'obs ^ 0.2. 
This would result in a rising light curve at early times, for all frequencies. On the other hand, 
for (a, b) = (2, 0), where the initial Lorentz factor is the same for all 9, the deceleration time, 
tdec; is very small everywhere, so that the initial rise of the light curve at t < tdec will be very 
hard to observe. The lack of an observation of a rising hght curve for afterglow observations 
starting from a few hours after the burst can already constrain the jet structure: either b < 2 
or Fq > 10^ are required. Future observations at much earlier times after the burst, as may 
be achieved with the forthcoming Swift mission, may provide much stronger constraints on 
the jet structure. 

It can also be seen from Figure 4 that for (a, b) = (2, 2), the value of the temporal decay 
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slope, a (defined by F,^ oc t"), before the jet break, is higher for larger viewing angles ^obs- 
This effect is very large for model 1, but not very significant for model 2 (see Figure 5). The 
value of a before the break in the light curve for model 1 is larger than the observed value 
for several well studied GRB afterglows, and the light curve power-law steepening, \Sa\, is 
smaller than the observed value for many GRBs. The light curves obtained using model '1' at 
early times are almost identical to the light curves from hydrodynamical simulations, and this 
suggests that jet structure described by (a,b) = (2,2) is not consistent with the observations 
for a number of GRBs. 

For (a, b) = (2, 0) the correlation between a before the jet break and 6obs is not seen, 
but instead there is a flattening of the light curve (i.e. an increase in a) just before the jet 
break, for viewing angles sufficiently large compared to the core angle of the jet, 9ohs ^ 3^c- 
This effect is more prominent in model 1, where it can also be seen to a lesser extent for 
(a, b) — (2, 2) as well, compared to model 2, where this effect is smaller and can be seen 
only for {a,b) = (2,0) (see Figure 5). Light curves calculated for the same jet profiles but 
with the jet dynamics given by a hydro-simulation (Kumar & Granot 2002) show that for 
(a, b) = (2, 0) a reasonably sharp jet break is obtained only for ^^obs ^ (2 — 3)6c, while the 
flattening in the light curve just before the jet break becomes prominent at ^obs ^ (5 — 7)6c. 
This leaves only a factor of ~ 3 in Oohs/^c ior which there is a sharp jet break with no 
flattening in the hght curve just before this break, as is seen in all afterglows with jet breaks. 
However, the inferred values of ^obs range at least one order of magnitude, ^obs ~ 2 — 20°. 
Therefore, (a, b) — (2, 0) is unhkely to be a universal jet profile, for all GRBs, since 9c 
is expected to vary between different afterglows. Even without requiring a universal jet 
profile, i.e. the same 6c for all GRB jets, it would be highly improbable if all viewing angles 
would by chance fall within the narrow range that gives a light curve qualitatively similar 
to observations. This can exclude (a, b) = (2, 0) for most GRBs, and imply 6 > (or even 
6 > 1). It is interesting to note that this provides a lower limit on 6, while the lack of 
detection of a rising light curve at early times may provide an upper limit on b. This may 
eventually pin down the value of b, or alternatively, call into question the simple jet structure 
interpretation discussed in this paper if the data is inconsistent with a reasonable value for 
b. 

Figure 6 shows the light curves for an external density that drops as {k = 2). This 
corresponds to a stellar wind of a massive star progenitor. In the pulsar wind bubble model of 
GRBs, the effective external density can have different profile which may be approximated as 
a power law with an index k that ranges between and 1 (Konigl & Granot 2002). Therefore 
k = and k = 2, are the extreme values of k that may be expected, and an intermediate 
value of k, resulting in an intermediate behavior of the light curves, is possible. 
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For k — 2 we see a, break in the light curve for a spherical flow {a — b — 0), at t ~ 10 yr, 
corresponding to the non-relativistic transition time, tNR- This can also be seen in Figure 
7, which shows the temporal decay slope, a, for the same light curves shown in Figure 6. 
Generally, £'1/(3-'=)^ so that for A; = 2 it is linear in E, and a more moderate total 

energy of, say, 10^^ erg would give tNR ~ 1 yr- We extend the observer time for which 
we show the light curves in order to illustrate that in the non-relativistic regime, all light 
curves become similar. One should keep in mind, that physically, one may expect an 
external density profile only up to some finite radius (for a stellar wind, this correspond to 
the radius of the wind termination shock) while at much larger radii we expect a roughly 
constant density medium. This would result in modifications to the light curves shown here, 
however, this is typically expected to effect the light curves only in the non-relativistic regime 

The deceleration time, tdec, is typically very small for A; = 2 since the density at small 
radii is very large, and the total swept up mass is only linear in radius. The jet break is much 
smoother and less sharp compared to a constant density environment {k — 0), as was found 
for a 'top hat' jet considered by Kumar & Panaitescu (2000). There is a relatively sharp 
break only for (a, 6) = (2,0), and even then it is hardly sharp enough to reproduce the jet 
breaks observed in GRB afterglows. This suggests A; < 2, or rather k < 1. It is interesting 
to note that this is consistent with the external density profile expected in the pulsar wind 
bubble model (Konigl & Granot 2002; Guetta & Granot 2002), while this is not consistent 
with the density profile expected in the coUapsar model, k — 2. 



4. Discussion 

We have calculated light curves from structured relativistic jets (jets whose Lorentz 
factor and energy per unit solid angle vary smoothly with angle from the jet axis) using two 
simple models for the jet dynamics, motivated by the hydrodynamical simulation of Kumar 
and Granot (2002), and as two limiting cases of energy redistribution in the lateral direction. 

The first model considers the energy per unit solid angle, e{9, t), to be time independent, 
so that each segment of the jet evolves independently of the other parts of the jet, as if it 
was part of a spherical fiow. This is roughly consistent with the results from hydrodynamical 
simulations of structured jets for which e varies slowly with 9, where we find the transverse 
velocity in the comoving frame to be small compared to the sound speed throughout much 
of the jet evolution, as long as the Lorentz factor, F, is of order a few or larger along the 
jet axis (Kumar & Granot 2002). This model is expected to be a good approximation for 
calculating the jet dynamics and the resulting afterglow light curves for the first few days 
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after a GRB. 

The second model considers the maximum possible re-distribution of energy in the 
transverse direction so as to reduce the lateral-gradient of e. The energy per unit solid angle, 
e{9,t), is taken to be proportional to the average over the initial distribution of e, e{9,tQ), 
within the area out to which a sound wave could have traveled from 9 since the initial time, 
to- This model is likely to be more accurate in describing the late time behavior of the jet 
and the light curves, when the transverse velocity becomes of order the sound speed. 

We have calculated light-curves for a number of different initial jet structures, which 

are taken to be power-law profiles at angles larger than some core angle, i.e. e oc 9^'^ and 
r oc 9~^ for 9 > 9c- We have considered three different jet profiles, namely {a,b) = (2,0), 
(2,2), (0,2) for both a homogeneous external medium as well an ambient density falling off 
as distance squared from the center of explosion (i.e. A; = or 2, where Pe^t o<: r~^). 

The quantitative differences in the light-curves calculated using our two models are 
typically of order unity for all of the jet structures we have considered. This gives us some 
confidence that lightcurves can be calculated with a reasonable accuracy even with a crude 
modeling of jet dynamics. There are, however, interesting qualitative differences in the light- 
curve properties depending on how we model jet dynamics. For instance, for (a, h) = (2, 2) 
and a homogeneous medium {k — 0), the temporal index of the light curve, a (defined by 
oc t"), before the jet break is larger for large viewing angles, ^obs- This effect is much 
more prominent for model 1, compared to model 2. For (a, b) — (2, 0) and k — 0, with model 
1 there is a pronounced flattening in the light curve just before the jet break for ^obs ^ 3 9c, 
while for model 2 this effect is much smaller. 

There is a significant difference between the two models in the late time light curves, as 
we expect: model 1 light curves continue a power-law decline at late times, whereas model 2 
light curves show a slight fiattening at late times resulting from energy redistribution - see 
Figures 4 and 5 for homogeneous external medium with jet structures (a, b) = (2, 2) & (2,0), 
and Figures 6 and 7 for the stratified external medium the particular case of (a, b) = (2, 0). 
Figures 6 and 7 also show that the break in the light curves for a stratified external medium 
is very gentle and the change to the power-law index a is small so long as the observer is not 
too far the axis: 9ohs/9c ^ 10. This result is consistent with the work of Kumar & Panaitescu 
(2000), based on the analysis of a 'top hat' jet model, that it is very difficult to see a break 
in light curves for jets in an r"^ medium. 

We can learn about the jet structure in GRBs from comparing our theoretical light 
curves to GRB afterglow observations. For instance, the light curves for (a, b) = (2, 0) and 
(2,2) behave very differently prior to the jet break. This can be used to discriminate between 
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these possibilities (Figures 4 & 5). For obvious reasons, the hght curves for the (0,2) case are 
very distinct and similar to a spherically symmetric explosion, and do not show a jet break. 
This jet structure can therefore be ruled out for all GRBs that show a jet break in their 
afterglow light curves. Clear breaks in the light curves can be seen for a stratified external 
medium density only when the observer is located far away from the jet axis (^obs ^ 10 ^c), 
however, in this case the GRB will be faint and is likely to be missed in flux limited triggers 
for detecting GRBs. 

The light curves we obtain for structured jets, in many cases show a different qualitative 
behavior compared to afterglow observations. For example, k = 2 docs not produce a 
sufficiently sharp jet break, while for k = only a ~ 2 produces jet breaks that resemble 
afterglow observations, and even then, it is very difficult to produce all current afterglow 
observations with a universal jet profile, that is, if all GRB jets have the same structure and 
differ only by our viewing angle, ^obs, w.r.t. the jet symmetry axis. Even if GRB jets do not 
have a universal structure, a comparison between the observed light curve and theoretical 
calculations can constrain the jet structure of each GRB separately. 

Future afterglow observations should enable us to determine the structure of relativistic 
jets in GRBs as well as the properties of the surrounding medium. This will enable us to 
determine if the energy release in GRBs is nearly constant or not, and whether the observed 
afterglow light curves and the jet break time arc determined by the observer angle w.r.t. a 
structured jet, or by the opening angle of a 'top hat' jet. 

This work was supported by the Institute for Advanced Study, funds for natural sciences 
(JG). 
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Fig. 1. — The evolution of the Lorentz factor r(^) — 1, for (a, 6) = (0,2), Fq = 10^ and 
6c = 0.02 rad. For a = 0, models 1 and 2 become the same, so that this figure applies to 
both models. 
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Fig. 2. — The hydrodynamic evolution of the energy per unit sohd angle e and the Lorentz 
factor r — 1, as a function of the angle 9 from the jet symmetry axis, according to our two 
simple and extreme models, for (a, 6) = (2,0). 
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Fig. 4. — Light curves for a constant density external medium {k = 0), for models 1 and 2 
(see text for details), in the optical (z/ = 5 x 10^^ Hz), for a jet core angle 9c = 0.02, viewing 
angles ^obs = 0.01, 0.03, 0.05, 0.1, 0.2, 0.3, 0.5, p = 2.5, = ee = 0.1, n=l cm'^, Tq = 10^, 
and eo was chosen so that the total energy of the jet would be 10^^ erg. A power law of 
is added in some of the panels, for comparison. 
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Fig. 5. — The temporal decay slope, a = din F^/ dint (i.e. F^, oc for the light curves 
shown in figure 4, that are with a constant density environment {k = 0). 
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Fig. 6. — The same as Figure 4, but for an external density oc r ^, corresponding to a stellar 
wind, with A = 5 x 10^^ gr cm~^. 
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Fig. 7. — The temporal decay slope, a, for the light curves shown in figure 6, that are for a 
stellar wind environment {k = 2). 



